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We present a now measurement of the Z/^" transverse momentum distribution in the range 
- 330 GeV, in proton-antiproton colhsions at = 1.96 TeV. The measurement uses 0.97 fb^^ of 
integrated luminosity recorded by the DO experiment and is the first using the Z/^" — >■ ^'^ ijT+X 
channel at this center-of-mass energy. This is also the first measurement of the .^/7* transverse 
momentum distribution that presents the result at the level of particles entering the detector, min- 
imizing dependence on theoretical models. As any momentum of the Z/^* in the plane transverse 
to the incoming beams must be balanced by some recoiling system, primarily the result of QCD 
radiation in the initial state, this variable is an excellent probe of the underlying process. Tests of 
the predictions of QCD calculations and current event generators show they have varied success in 
describing the data. Using this measurement as an input to theoretical predictions will allow for a 
better description of hadron collider data and hence it will increase experimental sensitivity to rare 
signals. 

PACS numbers: 12.38.Qk, 13.85.Qk 
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In the complex environment of a hadron collider, such 
as the Fermilab Tevatron Collider or the CERN Large 
Hadron CoUider, the — > e+e^ and — > ^+^^ 
processes are experimentally simple to identify and have 
little background. Further, reconstruction of the Zj^* 
kinematics provides an unambiguous, colorless probe of 
the underlying collision process. Momentum conserva- 
tion requires that any momentum component of the Z j^* 
in the direction transverse to the incoming hadron beams 
(j)t) must be balanced by a recoiling system (X), typ- 
ically the result of QCD radiation in the initial state. 
The Zj^* pt is therefore sensitive to the nature of this 
radiation across a wide momentum range, making it a 
compelling variable and an excellent testing ground for 
theoretical predictions. 

Several tools have been developed which give a predic- 
tion of the Z/ 7* Pt distribution, from fixed order pertur- 
bative QCD (pQCD) calculations valid at high px, such 
as MCFM FEWZ [3] and DYNNLO Q, to predictions 
based on gluon resummation valid at low pT [J] , such as 
RESBOS [5|. Various complete event generators are also 
available, including pythia herwig 0|, ALPGEN [1], 
and SHERPA [9|] , which cover both high and low px regions 
by interfacing tree-level matrix element calculations to 
a parton shower resummation model. Comparisons be- 
tween these generators show that they differ significantly 
in the predicted kinematics of boson and boson-|-jet pro- 
duction, and that these predictions have a strong depen- 
dence on various adjustable internal generator parame- 
ters Measurements of the Z/j* px and other kine- 
matic quantities in Z/j* production are therefore an es- 
sential input to improve these models, which are also used 
to predict the properties of rare signals like the Higgs bo- 
son and its main backgrounds: W-|-jets, ^/7*-|-jets and 
diboson production. Such improvements will result in 
increased experimental sensitivity to these rare signals. 

Previous measurements at the Tevatron have stud- 
ied the Z/7* Pt and rapidity distributions both inclu- 
sively (llHl^ l and in events with at least one jet [l5| . 
Other measurements have focused on the kinematics of 
the jets in Z/j* or W boson events 15-1^, of the an- 
gular correlations between the Z/ 7* and leading jet [lOl , 
and of the production of Z/j* or W boson in associ- 
ation with heavy flavor quarks (2ll - E^ . In this Letter, 
we describe a new measurement of the normalized inclu- 
sive Z/7*— >■ jX^iiT Pt distribution, the first such mea- 
surement using the dimuon channel in the Tevatron run 
beginning 2001 ("Run II"). The differential dimuon+X 
cross section is measured as a function of the dimuon pT 
(Pt), then normalized to the measured inclusive dimuon 
cross section, canceling many systematic uncertainties. 
The shape of the Z/j* pT distribution has previously 
been measured in Run II with the Z/j* — >■ e~^e~ chan- 
nel [l3| . using a comparable integrated luminosity. In 
that result, a resummation prediction was found to be 
consistent with the data in the pt < 30 GeV region. 



but pQCD predictions were found to be 25% below the 
data in the region pT > 30 GeV. Compared to that re- 
sult, the muon channel uses a statistically independent 
dataset, has a different detector acceptance, and differ- 
ent sources of systematic uncertainty; it therefore adds 
important information on the Z/^* pT distribution and 
any disagreements between the data and theory predic- 
tions. Due to the different response of the detector to 
electrons and muons, there is also different sensitivity to 
QED final state radiation (FSR) in the dielectron and 
dimuon systems 25 1. 

An important development in this analysis is the def- 
inition of the final observable: for the first time in a 
measurement of the Z/j* pT, the results are presented 
at the level of particles entering the detector. Previous 
measurements have applied theoretical factors to go from 
these particles to the (non-observable) Z/^ by correct- 
ing for any undetected FSR, and from the measured lep- 
ton acceptance to full Att coverage, correcting for unde- 
tected leptons. These factors rely upon models of FSR 
and the correlation between boson rapidity and pt- Here, 
we avoid such factors and present the data in terms of 
an observable: the pt of the dimuon system, for muons 
within the detector acceptance. This approach minimizes 
dependence on theoretical models, and the result can be 
used as an unbiased test of such models. This is also 
the same definition of the dimuon final state as previ- 
ous DO measurements of Z/7*(— )• fj,~^ fi~)+]et+X pro- 



duction [15|, |20[, and the relationship between p^ and 



the production of jets in the final state makes this mea- 
surement complementary to those results. 

The analysis uses a dataset of pp collisions at 
y/s = 1.96 TeV, corresponding to an integrated luminos- 
ity of 0.97 ± 0.06 fb-i [iel recorded by the DO detec- 
tor between April 2002 and February 2006. A fuU de- 
scription of the DO detector is available elsewhere 27| . 
and only a brief description of the components most rel- 
evant for this analysis is given here. The pp interaction 
region is surrounded by two tracking detectors: a sili- 
con microstrip tracker and a scintillating fiber tracker, 
both housed inside a solenoidal magnet providing a field 
of approximately 2 T. These trackers provide a momen- 
tum measurement for charged particles and are used to 
reconstruct the primary interaction point in each colli- 
sion. Outside the solenoid is a liquid-argon and uranium 
calorimeter which is split into three sections: a central 
section covering {rjl < 1.1 [28| and two forward sections 
covering 1.4 < |?7| < 4.2. Outside the calorimeter there 
are three layers of muon detectors, made of a combina- 
tion of scintillation counters and drift chambers covering 
1 77 1 < 1, and scintillation counters and drift tubes extend- 
ing the coverage to \r]\ < 2. A 1.8 T iron toroidal magnet 
is located between the first and second layer, providing 
an independent momentum measurement for muons. 

Events used in this analysis are selected by at least 
one of a suite of single-muon triggers. These triggers 
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used a fast readout from muon system, or a combination 
of the fiber tracker and muon system, to identify muon 
candidates. Then information from the full tracking and 
muon systems is incorporated to provide further rejec- 
tion. Additional requirements are then applied to the 
events selected by the trigger to obtain a sample of Zj^* 
candidates. Using the full information from the muon 
detectors and the tracking system, two muons of oppo- 
site charge and pt> 15 GeV are required, with a dimuon 
mass in the range 65 < < 115 GeV. To reject cosmic 
rays and poorly reconstructed muons, the muon tracks 
are required to be consistent with the reconstructed pri- 
mary interaction point both along the beam direction and 
in the transverse plane. The two muon tracks are also re- 
quired not to be colinear, and to be consistent with the 

bunch crossing time using timing information from the 
muon system scintillators. Further selections are applied 
to limit the measurement to regions with high detection 
efficiency: the muons are required to have \'r]\ < 1.7, and 
the primary vertex must lie within 50 cm of the center of 
the detector in the coordinate along the direction of the 
beam. In total, 59,336 dimuon candidate events pass all 
selection requirements. 

The main background in this analysis is dijet produc- 
tion with two semi-leptonic decays, or H^-|-jets produc- 
tion in which one muon comes from the W and the other 
from a semi-leptonic decay in a jet. These events are 
reduced to a negligible level by two isolation require- 
ments. First, we reject overlaps between muons and 
jets with pt> 15 GeV, by requiring angular separation 
-y/ (A0)^ -I- (At;)^ > 0.5, where </) is the azimuthal angle. 
Then, we require the product of the isolation variables for 
the two muons to be < 0.05, where each isolation variable 
is calculated by taking the sum of track and calorime- 
ter energy in a cone around each muon (excluding the 
muon track and calorimeter energy associated with the 
muon itself), and dividing by the muon pT- 

The remaining contribution from these backgrounds is 
estimated from data by studying the product of the iso- 
lation variables for muons failing isolation requirements. 
Extrapolating into the selected region shows this back- 
ground to be < 0.5% of the final sample. The remain- 
ing backgrounds (from tt, WW, WZ, and Z/j* t+t^ 
production), as well as the Z/j* — > fi^^^+X sig- 
nal, are modeled with pythia v6.409. A separate 
Z/'y* i.i^fi~+X sample is generated using ALPGEN 
v2.11 with PYTHIA v6.409 for parton showering. All 
simulated signal and background samples are normal- 
ized to higher order theoretical predictions 2^ 3] and 
passed through a GEANT [31[ simulation of the DO detec- 
tor. The total background from all sources is found to 
be below 2% everywhere and less than 1% in the region 

< 50 GeV. The estimated background contribution is 
subtracted from data and a 10% systematic uncertainty 
is assigned to each background normalization to cover all 
sources of uncertainty. 



To extract the shape of the distribution, the mea- 
sured dimuon candidates must be corrected for detector 
resolution and efficiency, both of which are derived di- 
rectly from data. The detector resolution is extracted 
from the shape and position of the Z/j* resonance peak 
in dimuon data, which is dominated by detector resolu- 
tion rather than the natural width of the Z boson. The 
resolution is well described by a double- Gaussian func- 
tion form in l/px, with the majority (98%) of muons 
having a 1 /pt resolution of approximately 0.0018 GeV~^ , 
and the remaining 2% (chosen at random) a 1/pt reso- 
lution of approximately 0.012 GeV^^. 

The detector efficiencies are derived using the "tag and 
probe" method on dimuon candidate pairs. The "tag" 
muon is selected, and must pass all selection require- 
ments, which may be adjusted as needed to remove back- 
grounds. The "probe" muon is then selected with one ex- 
plicit reconstruction requirement removed; the fraction of 
probe muons which also meet this requirement gives an 
measurement of the efficiency. In this way, the efficiency 
of reconstruction, trigger and isolation requirements are 
measured individually, and parameterized in terms of the 
geometry of the detector. The method is repeated on 
simulated events, where results typically agree with data 
to within 3%, and factors are applied to the simulated 
events to correct for any such discrepancies. However, 
the muon trigger is not simulated; instead the trigger ef- 
ficiency measured in data, with an average efficiency of 
approximately 88%, is applied on an event-by-event basis 
to the simulated events. 

The binning used for the data is selected based on 
a combination of detector resolution and data statistics 
considerations. Detector effects on the resulting pf dis- 
tribution are assessed by comparing the p^ defined in 
terms of particles entering the detector to the pf recon- 
structed by the detector. First, we define the dimuon sys- 
tem at the particle level in a way that can be implemented 
in any simulation. We consider all particles with lifetimes 
> 10 ps to have reached the detector. From this list of 
particles, all muons with px > 15 GeV and jryj < 1.7 are 
selected (regardless of their source in the generator event 
record), matching the detector acceptance. By construc- 
tion, the muons are considered after QED FSR, as would 
be measured in the tracking detector. Then, all possible 
opposite charge muon pairs are formed, and any which lie 
within the required mass range of 65-115 GeV are kept. 
In the rare cases (< 0.5%) of events with more than one 
selected pair, the pair with mass closest to the Z bo- 
son mass is used. The requirement of non-colinearity of 
muon tracks applied at detector level is found to reject 
less than 0.2% of particle level candidates, so is not ap- 
plied at the particle level. Similarly, the muon isolation 
and vertex requirements are treated as a detector level 
selection which is corrected for, and are not implemented 
at the particle level. 

We next correct the measured p^ in data to the parti- 
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cle level, using the alpgen+pythia Z/7* 
sample. There are three possible scenarios for any given 
event, resulting in a three step process. In the first sce- 
nario, a dimuon pair may pass all detector level selec- 
tions, but fail one or more particle level selections. This 
class of events is dominated by migrations into the se- 
lected mass or muon {rjl regions due to detector resolu- 
tion. In the simulated events, this class makes up approx- 
imately 2% of the final sample with negligible dependence 
onp^, and this predicted contribution is subtracted from 
the measured data. In the second scenario, the dimuon 
pair passes both the particle and detector level selections. 
For these simulated events, the particle level is plot- 
ted against the detector level to assess the impact of 
detector resolution (see Fig[T]). The data distribution af- 
ter the subtraction described in the first step is then cor- 
rected using a regularized inversion of this resolution ma- 
trix 3^. The regularization imposes the condition that 
second derivatives be small, which produces a smooth 
distribution; this smoothing is accounted for when deriv- 
ing the uncertainties. In the third scenario, the particle 
level dimuon pair may pass selections, but the detector 
level pair fail selections. This effect is also assessed using 
simulated events and is dominated by inefficiency in the 
trigger or reconstruction and by gaps in detector cover- 
age within the muon acceptance. The data distribution 
resulting from the second step is corrected for these in- 
efficiencies, parametrized as a function of (particle level) 
, giving the Z/^*+X differential cross section. 
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FIG. 1: The resolution matrix for from simulated events, 
used in correcting the measured data to particle level. The 
box area represents the number of events in a given particle 
level bin that are in a given bin of measured (detector 
level) py. 

Ensemble testing is used to determine the uncertainty 



on the differential cross section and any biases in the 
correction process. The PYTHIA Z/j* — !• fi'^fj,~+X sam- 
ple is used to build pseudo-datasets, after first being 
adjusted to describe the data by applying the ratio of 
the fully corrected data to the particle level in 
this PYTHIA sample. Three hundred pseudo-datasets are 
drawn, with events chosen at random for each pseudo- 
dataset with a probability set so the average pseudo- 
dataset size matches the measured dataset. Each pseudo- 
dataset is then treated exactly as data, and the detector 
level distribution corrected using the three step process 
described above. The resulting corrected distribution is 
compared to the true particle level in that pseudo- 
dataset, and the fractional difference (r,) is calculated 
for each p^ bin. This process is repeated for all 300 
pseudo-datasets. In a given p^ bin, the 300 form a 
Gaussian distribution and any shift of the mean of this 
distribution away from zero indicates a bias in the cor- 
rection process. The data are corrected for such biases, 
which are all at or below the 1% level, and the uncer- 
tainty on the Gaussian mean is assigned as a system- 
atic uncertainty. The RMS of the distribution re- 
sults from comparing the corrected distribution for each 
pseudo-dataset, which is smoothed by the regularization 
in the matrix inversion step, to the true particle level for 
that pseudo-dataset, which contains statistical fluctua- 
tions. This RMS is therefore assigned to the data points 
as the statistical uncertainty. It is comparable to, but al- 
ways larger than, l/V-^data, where A'data is the number 
of detector level data events in a given bin. 

Finally, further systematic uncertainties are assessed. 
Several sources are considered, including the muon mo- 
mentum resolution and the absolute momentum scale, 
and various detector efficiencies for muons. Each source 
is shifted individually up and down by one standard de- 
viation in the ALPGEN-I-pythia sample used to derive 
the correction factors, and left unchanged in the pseudo- 
datasets. The full ensemble testing process is repeated. 
Statistical uncertainties are unaffected by these tests, but 
any change in the correction bias is assigned as a system- 
atic uncertainty. All systematic uncertainties are com- 
bined in quadrature, with the muon 1/pt resolution be- 
ing the largest source at low p^. At high p^, detector 
efficiency effects are the largest systematic source, but 
the statistical uncertainties dominate. 

Finally the differential cross section is normalized to 
the total dimuon production cross section (with the same 
muon 1 77 1, pt and dimuon mass requirements) measured 
in the data, determined by integrating over all p^. The 
dominant systematic uncertainties on the total cross sec- 
tion arise from the absolute determinations of the lumi- 
nosity and muon trigger efficiency. Neither of these has a 
dependence on p^ , and therefore they do not contribute 
to the uncertainty on the normalized distribution. 

Table [J lists the normalized differential cross section, 
(l/cr) X (dcr/dp^), together with the statistical and sys- 
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tematic uncertainties. We note that, due to the smooth- 
ing introduced by the regularization condition imposed 
during the second step of the data corrections, statistical 
fluctuations in the measured cross section in each bin 
have been suppressed; however the statistical uncertain- 
ties still accurately reflect the possible spread in each bin 
which could be caused by such fluctuations. As a result, 
care must be taken when using the data in any fits as this 
suppression of fluctuations may lead to an artificially low 
for any fit which describes the central values of the 
data well. Table |T] also lists four multiplicative correction 
factors for each bin, which can be applied to compare this 
result to previous measurements: the factor labeled 
corrects for the effect of the muon pt> 15 GeV require- 
ment; the factor labeled FSR corrects for QED FSR; the 
factor labeled A then corrects from the measured Icpton 
acceptance to full An acceptance; and finally the factor 
labeled M corrects from the measured mass window to 
the larger mass window used in the DO electron channel 
measurement (40-200 GeV). Applying only the 
factor results in the same dimuon definition as previous 
^/7*+jets measurements [l^ E^; unlike p^, the vari- 
ables studied in these previous measurements had mini- 
mal dependence on the muon px requirement, so a cor- 
rection was applied by default. All factors are derived 
using RESBOS interfaced to PHOTOS [s^], as described in 
the following text, and we provide only the central val- 
ues without assessing possible systematic uncertainties. 
However, deriving the same factors from the different the- 
oretical calculations described in the following text indi- 
cates that model dependence limits the accuracy of these 
factors, particularly A for p^ > 20 GeV, to the level of 
a few percent. Applying all factors to the data allows a 
comparison to the DO electron channel measurement, as 
shown in Fig. [2l Within the limitations of this compari- 
son, the agreement is reasonable. For direct comparisons 
with theory, these correction factors are not applied to 
the dimuon data. 

To compare to the data, predictions for the p^ dis- 
tribution are obtained from several theoretical calcula- 
tions. Predictions from pQCD are obtained with mcfm, 
by evaluating both the differential distribution and to- 
tal cross section at either leading order (LO) or next-to- 
leading order (NLO): 
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where the first term is of order 0{a^) at LO and 0{al) 
at NLO, while the second term is 0{al) at LO and 
0{ag) at NLO. This approach differs from the treat- 
ment of the pQCD calculation in the DO electron chan- 
nel measurement. There, both the total cross section 
and differential distribution were calculated to the same 
power of the strong coupling constant, 0{a'g), yield- 
ing a NNLO total Z/j* cross section (and was labeled 
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FIG. 2: Measurements of the normalized differential cross 
section in bins of for the dielectron [T^] and dimuon chan- 
nels. Both results are shown with combined statistical and 
systematic uncertainties. 



"NNLO"), but a NLO differential distribution. The pre- 
scription used here, calculating both the total and differ- 
ential cross section to the same number of contributing 
terms in the perturbative expansion, results in a reduced 
scale dependence and improved convergence of the per- 
turbative series 



37j . The total cross section is evalu- 



ated using the inclusive pp ^ Z/j* + X process at LO 
and NLO, and the differential distribution evaluated us- 
ing the pp — )■ Z/j*+ jet+X process again at LO and 
NLO, with no limit on the jet rapidity but requiring jet 
Pt > 2.5 GeV to remove the divergence as pt — t- 0. The 
same requirements are placed on the muons as for the 
data analysis, and the differential Z/^* pT distribution 
close to the jet px cutoff is excluded. The MSTW2008 
LO and NLO PDFs [s^ are used throughout in calcu- 
lating the LO and NLO processes respectively. In all 
cases, renormalization and factorization scales are set to 
the sum in quadrature of the mass and pt of the Z/"f* in 
each event, and the dependence on this choice is assessed 
by varying both scales simultaneously up and down by a 
factor of 2, both for the differential distribution and the 
inclusive Z/"f* cross section used in normalization. PDF 
uncertainties are assessed using the MSTW2008 68% er- 
ror sets, again taking into account the effect on the dif- 
ferential distribution and the inclusive Z/j* cross section 
used in normalization. These are found to be approxi- 
mately a factor of two smaller than the scale uncertainties 
at NLO, and negligible compared to the scale uncertain- 
ties at LO. The prediction from MCFM must then be cor- 
rected for the effects of QED FSR from the muons. These 
corrections are derived from the RESBOS-f PHOTOS sample 
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described below, by comparing the distribution obtained 
by constructing the dimuon state using the muons before 
and after QED FSR, and are the inverse of the per-bin 
FSR corrections hsted in Table HI We note that apply- 
ing the pQCD prescription from the DO electron channel 
measurement here would change the LO prediction to 
LO/NLO, and the NLO prediction to NLO/NNLO. Dif- 
fering only in the total cross section used to normalize, 
this change would lower the current predictions by 28% 
and 4% respectively across all Zj^* pr, and increase the 
scale uncertainties by ~ 10% (a < 1% absolute increase 
in the uncertainty). 

A prediction for (1/cr) x (da/dp^) is also obtained 
from RESBOS, using the CTEQ6.6 parton distribution 
functions (PDF) [3J|. At low , RESBOS performs a 
next-to-next-to-leading logarithm (NNLL) resummation 
calculation, using the BLNY parametrization [sHj, with 
the default settings taken here. At higher p^ {p^ > 
20 GeV), RESBOS transitions to an 0{as) pQCD calcula- 
tion. Across the entire p^ range, k-factors are applied in 
bins of Z/^* pt and rapidity to normalize to an 0{al) 
pQCD calculation [36|. The renormalization and factor- 
ization scales for the pQCD calculation are set to the 
mass of the Z/j* in each event. RESBOS is interfaced to 
PHOTOS for the simulation of QED FSR. 

Predictions for (1 /a) x (dcr/dpf ) are also obtained from 
four event generators. Previous measurements [isl [Til . 
[ioj indicate that the best description of boson+jcts fi- 
nal states is currently provided by LO 2 ^ N matrix 
element calculations with matched parton showers, as 
implemented in SHERPA and ALPGEN, so we focus on 
these. We use the same PDF set for all event genera- 
tors: CTEQ6L1 [il]. First, a sample of events is gener- 
ated with SHERPA, which uses the COMIX matrix element 
generator interfaced to a Catani-Seymour subtrac- 
tion based parton shower [41| . Here, up to three partons 
are included in the matrix element calculation, and the 
threshold for matching to the parton shower is set to the 
default value of 30 GeV. In SHERPA, the scales are deter- 



mined dynamically during the matching process 42|. A 
sample of events is then generated with ALPGEN, again 
with up to three partons in the matrix element calcula- 
tion. The factorization scale is set to the sum in quadra- 
ture of the mass and pt of the Z/j* in each event, and 
the renormalization scale set according to the CKKW 
prescription [43| . Parton jets from the matrix element 
calculation are required to have pt> 13 GeV, Ai?(jet, 
jet) > 0.4, and are limited to |?7| < 2.5. These events 
are hadronized in three ways: first, using herwig (using 
an angular ordered parton shower) with jimmy [44| for 
multiple parton interactions, then using PYTHIA with un- 
derlying event tune D6 45| (using the virtuality-ordered 
shower), and finally using PYTHIA with tune Perugia 
6 [i^ (using the pT-ordered shower 47] ) . This results in 
three different ALPGEN predictions. In each case the de- 
fault matching procedure is applied after hadronization. 



requiring a Ai?(jet, jet) < 0.4 match between parton jets 
and particle jets with pT > 18 GeV. To determine the 
impact of the matching to the ALPGEN matrix elements 
calculation, herwig and pythia are also tested directly 
in the same configurations described above: herwig with 
JIMMY for multiple parton interactions, pythia with 
tune D6, and pythia with the Perugia 6 tune. In these 
configurations, all final state partons are generated by 
the parton shower. The renormalization and factoriza- 
tion scales for the hard scatter are set to the mass of 
the Z/j* in each event, and are determined dynamically 
for the initial and final state showers. For the RESBOS, 
SHERPA, ALPGEN, PYTHIA and HERWIG Calculations, the 
particle level quantities are extracted as defined earlier 
and each differential cross section prediction is normal- 
ized to the prediction of the dimuon cross section (with 
the same muon px, \ri\ and dimuon mass requirements) 
from that same model. 

The normalized differential cross section is presented 
in Fig. [3l The data points are placed at the bin av- 
erage, defined as the point where the differential cross 
section within the bin, taken from pythia reweighted to 
match the shape in data, is equal to the measured value 
in the bin [i^. For clarity, only the predictions of NLO 
pQCD and pythia Perugia 6 are shown with the data in 
Fig. ISJa). In the other parts of Fig. [Sj ratios are shown. 
To avoid repeating the data uncertainties and statistical 
fiuctuations several times, we do not use data as the de- 
nominator in these ratios. Instead, we choose pythia 
Perugia 6, as this provides the best overall description of 
the data, simplifying the determination of trends in other 
theoretical predictions relative to the data. As an exam- 
ple of the scale uncertainty in an event generator, two 
further pythia Perugia 6 samples are generated, with 
the scale for initial state QCD radiation varied up and 
down by a factor of 2. The effect of this change is shown 
as a shaded band around unity, and shifts the distribu- 
tion in opposite directions at low and high p^ , with the 
transition point at approximately 6 GeV. Further, even 
though PYTHIA is based on LO matrix elements, the scale 
uncertainty obtained is comparable to that on the NLO 
pQCD calculation, suggesting a cancellation of some of 
this scale variation in pythia through a detailed balance 
between the matrix element for Z/^* production and the 
Sudakov form factors from the parton shower. However, 
this small uncertainty is somewhat deceptive, as pythia 
does not include a full NLO calculation. 

Comparisons to the data indicate two regions: p^ < 
30 GeV, where the resummation calculation provides a 
good description of the data, and > 30 GeV, where 
the fixed order calculation provides the best description. 
In this higher p^ region, the NLO pQCD calculation is a 
significant improvement in uncertainty over LO, however 
an overall normalization difference relative to the data is 
observed. For 30 < < 100 GeV, this difference is be- 
tween 1-2 standard deviations of the combined data and 
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FIG. 3: The normalized differential cross section in bins of for Z/^* (— >■ fifi) + X events. The data are shown with statistical 
uncertainties (horizontal bar) and combined statistical and systematic uncertainties (full bar). The distribution is shown in (a) 
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parton shower generators in (d). All ratios in (b), (c), and (d) are shown relative to pythia Perugia 6. 



theory uncertainties, with the theory scale uncertainty 
dominating and the choice of a lower scale bringing the 
pQCD calculation into better agreement with the data. 
For > 100 GeV, the data statistical uncertainty dom- 
inates, and the theory remains below the data but is con- 
sistent within this uncertainty. This disagreement with 
pQCD predictions is in the same direction as observed in 
the DO electron channel measurement, but significantly 
smaller. However, the detector acceptance for the elec- 
tron channel was larger than for the muon channel, and 
attempting to extrapolate between these acceptances re- 
vealed a dependence on the theoretical models used, com- 
plicating direct comparisons of the two results. These 
two measurements in fact provide different information 
on the distribution in different Z/j* rapidity ranges, 
and future measurements which further probe the cor- 
relations between the Z/"f* px and rapidity are clearly 
of interest and may further illuminate the disagreements 
seen when comparing pQCD to data. 

Of the event generators, pythia Perugia 6 provides 
the best description of the data over the full range. 



and we note that the DO electron channel measure- 
ment was used as an input in deriving this tune. All 
other event generators agree within the combined the- 
ory and data uncertainties, except the pythia D6 and 
ALPGEN-I-PYTHIA D6 predictions. Interfacing pythia 
and HERWIG with ALPGEN clearly affects the region domi- 
nated by the ALPGEN matrix elements, though the agree- 
ment with data is equally good. While the SHERPA pre- 
diction agrees with the data within uncertainties, it gen- 
erally follows the shape of the pythia Perugia 6 predic- 
tion with a higher scale choice suggesting that, as for the 
pQCD calculation, a lower scale choice in SHERPA may 
yield an even better description of the data. 

In summary, we have presented a new measurement of 
the normalized Z/7*(— >■ cross section, differential 

in the dimuon pT- This is the first such measurement at 
the level of particles entering the detector, allowing un- 
biased tests of theoretical predictions. The measurement 
was made using a sample corresponding to 0.97 fb~^ of 
integrated luminosity recorded by the DO experiment in 
pp collisions at y^s = 1.96 TeV. The current best pre- 
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dictions for vector boson production at hadron colliders 
were tested, and these predictions have varied success in 
describing the data. In particular, the disagreement with 
pQCD seen in the electron channel measurement at high 
Z/7* pt is smaller within the kinematic acceptance of 
this measurement, and the use of a lower scale within 
the calculation further reduces the disagreement. An ac- 
curate description of both the low and high regions 
is also essential in predicting the production rates and 
kinematics of jets in association with the Z/'y*, and this 
result is an important input for the tuning of theoretical 
predictions. Improving the modeling of this process will 
lead to increased sensitivity of searches for rare and new 
physics. 
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TABLE I: The measured normalized cross section in bins of dimuon pr (pt) for Z/y* jX^ ix~ +X events. For each bin, 

wc present the bin boundaries, the bin average ({Pt), defined in the text), the normahzed difi^erential cross section, and the 
uncertainties. Uncertainties are split into statistical, systematic with no correlations between bins ("uncorr. unc"), and sources 
of systematic uncertainty that are correlated across all bins. The sources of correlated uncertainty are: 1) muon l/pr resolution, 
2) muon detection efficiency, 3) efficiency of all other selections. Factors to correct the result for the muon pT requirement 
(Pt)' QED FSR, muon acceptance (A), and to the larger mass window of 40-200 GeV used in the DO electron channel analysis 
(M) are also provided (defined in the text). 

p| (^^f) 100 X Stat. Uncorr. Source 1 Source 2 Source 3 pt^ FSR A Jd' 



(GeV) 


(GeV) 


(1/GeV) 


unc. (%) 


unc. (%) 


(%) 


(%) 


(%) 










0.0 


4.0 


1.5 


6.13 


±1.2 


±0.1 


T7.2 


±1.4 


±1.1 


0.993 


1.06 


1.05 


1.02 


4.0 - 


8.0 


5.9 


6.27 


±1.0 


±0.1 


Tl.5 


±0.6 


±1.2 


0.993 


1.01 


1.03 


1.01 


8.0- 


12. 


9.9 


3.84 


±1.2 


±0.1 


±4.3 


±1.1 


±1.3 


0.995 


0.97 


1.01 


1.00 


12. - 


16. 


14. 


2.50 


±1.3 


±0.1 


±5.9 


±0.7 


±1.1 


0.998 


0.95 


0.99 


0.99 


16. - 


23. 


19. 


1.43 


±1.3 


±0.1 


±5.4 


±0.6 


±1.1 


1.003 


0.98 


0.97 


0.98 


23. 


32. 


27. 


0.704 


±1.5 


±0.1 


±4.2 


±0.8 


±1.1 


1.012 


0.99 


0.94 


0.97 


32. 


45. 


38. 


0.332 


±1.9 


±0.1 


±2.7 


±1.1 


±1.1 


1.026 


1.00 


0.90 


0.96 


45. - 


60. 


52. 


0.150 


±2.9 


±0.2 


±1.4 


±1.6 


±1.1 


1.042 


1.01 


0.81 


0.95 


60. - 


80. 


69. 


0.0611 


±3.8 


±0.3 


±0.2 


±2.3 


±1.1 


1.059 


1.02 


0.75 


0.95 


80. - 


110. 


93. 


0.0203 


±6.0 


±0.6 


=F0.6 


±2.7 


±1.6 


1.081 


1.03 


0.67 


0.94 


110. 


150. 


130. 


0.00530 


±11.0 


±1.1 


=Fl.2 


±3.3 


±2.2 


1.121 


1.03 


0.60 


0.94 


150. - 


200. 


170. 


0.00116 


±20.7 


±1.4 


Tl.7 


±3.7 


±2.4 


1.143 


1.04 


0.55 


0.95 


200. - 


330. 


240. 


0.000123 


±42.0 


±2.7 


T2.1 


±4.1 


±2.6 


1.131 


1.05 


0.52 


0.95 



